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Abstract 

The quadrupole deformation properties of the ground and low-lying excited states 
of the even-even Magnesium isotopes with N ranging from 8 to 28 have been studied 
in the framework of the angular momentum projected generator coordinate method 
with the Gogny force. It is shown that the N=8 neutron magic number is preserved 
(in a dynamical sense) in 20 Mg leading to a spherical ground state. For the magic 
numbers N=20 and N=28 this is not the case and prolate deformed ground states are 
obtained. The method yields values of the two neutron separation energies which 
are in much better agreement with experiment than those obtained at the mean 
field level. It is also obtained that 40 Mg is at the neutron dripline. Concerning 
the results for the excitation energies of the 2 + excited states and their transition 
probabilities to the ground state we observe a good agreement with the available 
experimental data. On the theoretical side, we also present a detailed justification 
of the prescription used for the density dependent part of the interaction in our 
beyond-mean- field calculations. 



1 Introduction 



Thanks to the new breed of Radioactive Isotope Beam (RIB) facilities set up in 
the last years all over the world the study of exotic systems far away from the 
valley of stability has become one of the main topics in today's nuclear physics. 
On the neutron rich side, because neutrons do not carry electric charge, the 
dripline can extend far away from the valley of stability. Many neutron rich 
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regions have been the subject of detailed studies in recent years and just to 
give an example one should mention the regions around the magic neutron 
numbers N=20 and N=28. These exotic systems show many new interesting 
features which deserve a careful study [1-10]. Among them, the weakening of 
the neutron magic numbers N=20 and N=28 in some nuclei is probably the 
most challenging to study. 

It was in the neutron rich region N 20 where the breaking of a shell clo- 
sure was first detected [11] and related with a shape transition to prolate 
deformed shapes [12]. Ground state deformations in this region easily explain 
the anomalous isotope shift observed in 31 Na [13] as well as the decrease in 
the two neutron separation energy S2N in 31 > 33 > 35 Na and 30 Ne [11,13]. A large 
deformation has also been inferred in the N=20 nucleus 32 Mg through the 
measurements of B(E2,0i — > 2f) transition probabilities [14]. In addition, 
the low excitation energy of the first 2+ excited state [15,16] in 32 Mg and the 
ratio E(4 + )/E(2 + ) = 2.6 are consistent with the expectations for a rotational 
state. Recently, evidence for a rotational band in 31 Na has also been obtained 
[17]. Although indirect, another evidence for the deformed nature of the 32 Mg 
ground state comes from the strong ground state deformation of 34 Mg inferred 
from the large value of the B(E2, Of — > 2+) transition probability measured 
[18] in this nucleus and the E(A + )/E(2 + ) ratio of 3.18 [19,20] which is very 
close to the rotational limit. Further evidence for deformed ground states in 
the neutron rich Mg isotopes have also been obtained experimentally [21-23]. 
In addition, the quadrupole collectivity in 32_38 Si was experimentally stud- 
ied in [24]. Finally, an extensive account of older experimental results can be 
found in Refs. [25,26]. 

In the framework of the mean field approximation the 32 Mg ground state has 
been found to be spherical (see for example [27-32]) but it has also been 
argued that dynamical correlations may play an important role in this and 
other nuclei of the region. In particular, the results presented in [12,33-39] gave 
initial hints concerning the importance of the restoration of symmetries and 
configuration mixing in some of the exotic nuclei considered. From the Shell 
Model point of view, the large deformations around N=20 have been explained 
by invoking excitations of pairs of neutrons across the N=20 shell gap (2 hu 
configurations). In a small set of nuclei (the island of inversion) the coexisting 
2 hu intruder states actually fall below the normal hu states and become 
the ground state configurations. In this way Shell Model calculations with 
different levels of complexity [40-46] have been able to explain the increased 
quadrupole collectivity around N=20. 

The existence of anomalies in shell closures, as already found in nuclei with N 
~ 20, is one of the main features for the experimental and theoretical studies 
carried out up to now in the neutron rich nuclei with N« 28. The /^-unstable 
nuclei around 44 S have attracted particular interest since these nuclei play an 
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important role in the nucleosynthesis process. The experimental study of the 
/9-decay properties of 44 S (N=28) indicated [47,48] that this nucleus is de- 
formed. This has been confirmed in subsequent intermediate energy Coulomb 
excitation studies [49-51]. More recent experiments on nuclei of this region 
[52,53] suggest that 44 S is indeed a deformed nucleus but with strong shape 
coexistence making more challenging its theoretical description. 

The region N as 28 has been studied using the Skyrme Hartree-Fock (HF) 
model and the relativistic mean field approximation (RMF) [54-58]. These 
studies also showed the onset of deformation for some nuclei with N w 28. 
Surprisingly, in none of these mean field calculations the rotational energy 
correction has been considered in spite that it is known to be a key ingredient 
for a proper description of the energy landscape as a function of the quadrupole 
deformation. However, in the mean field calculations of [34] with several pa- 
rameterizations of the Skyrme interaction, in the calculation reported in [38] 
in the framework of the Bohr Hamiltonian method with the Gogny interaction 
or in our recent study [59] of the phenomenology of quadrupole collectivity 
around N=28 the rotational energy correction was explicitly taken into account 
either in an approximate way or exactly. The conclusions extracted from those 
calculations concerning the erosion of the shell closure at N=28 do not differ 
qualitatively (but they do in the quantitative side) from the ones extracted 
from the calculations without the rotational energy correction. This indicates 
that the N=28 shell closure can be more easily broken than the N=20 one (see 
also [60]). From the Shell Model point of view an erosion of the N=28 shell 
closure in the sulphur isotopes has also been found in calculations with the 
full sd shell for protons and the full pf shell for neutrons (see, for example 
[61]). 

The purpose of this paper is to perform a systematic study, in the nuclei 
20_40 Mg, of the quadrupole deformation properties of the ground and excited 
states. This isotopic chain contains the three spherical magic numbers N=8, 
N=20 and N=28 and therefore it is a good testing ground to examine both 
the systematic of deformation and the possible erosion of spherical shell clo- 
sures. Besides, most of the considered nuclei are examples where the mean 
field energy landscape as a function of quadrupole deformation shows pro- 
late and oblate minima which are practically degenerate in energy. Therefore, 
the correlation energies associated with the restoration of broken symmetries 
(mainly the rotational symmetry) and the collective quadrupole motion have 
to be considered at the same time. In the present study both effects are taken 
into account in the framework of the Angular Momentum Projected Generator 
Coordinate Method [37,39,59,62]. The main reason to perform an exact angu- 
lar momentum projection is that the usual approximations, like the Gaussian 
Overlap Approximation (GOA), used to compute the rotational energy cor- 
rection as well as the B(E2) transition probabilities are not expected to work 
well [36,63] at least for some of the nuclei considered in this study. 
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We have used in our calculations the Gogny interaction [64] with the parame- 
terization D1S [65]. The use of the Gogny interaction in this systematic study 
is justified not only by the fact that this interaction provides reasonable re- 
sults for many nuclear properties all over the nuclear chart, but also by the 
good description of the phenomenology of quadrupole collectivity in the re- 
gions N 20 and N ps 28 obtained recently in the same framework as the 
one used in this study [36,37,39,59] as well as in the context of the Bohr col- 
lective Hamiltonian [33,38]. As the results presented in this paper will show, 
this force also provides reasonable results for the systematic description of the 
phenomenology of quadrupole collectivity in the chain of Magnesium isotopes. 

The paper is organized as follows: in Section 2, we describe the theoretical 
formalisms used in the present paper. These formalisms are the Angular Mo- 
mentum Projection (AMP) and the Angular Momentum Projected Generator 
Coordinate Method (AMPGCM) [36,37,39,59,62] with the axially symmet- 
ric quadrupole moment as generating coordinate. The angular momentum 
projection technique is presented in section 2.1 where also special attention 
is paid to the specific problems appearing in the case of a density depen- 
dent interaction like the Gogny force. In section 2.2 we present the simplified 
expressions corresponding to the results of 2.1 in the case of axially sym- 
metric intrinsic wave functions. Configuration mixing of angular momentum 
projected Hartree-Fock-Bogoliubov states is presented in section 2.3. There 
we also present the expressions to compute both transition probabilities and 
spectroscopic quadrupole moments (for more details the reader is referred to 
Appendix A). The results of our calculations for several Magnesium isotopes 
are discussed in Section 3. In section 3.1 we present the results for the under- 
lying mean field studies. In section 3.2 the topological changes introduced in 
the mean field potential energy surfaces due to the exact restoration of the 
rotational symmetry are presented and then the results of the angular momen- 
tum projected configuration mixing calculations are discussed and compared 
with the available experimental data and other theoretical calculations. Fi- 
nally Section 4 is devoted to some concluding remarks. 



2 Theoretical framework: Angular momentum projection and con- 
figuration mixing 

As the results of the next sections will show, the mean field approximation 
is just an starting point and additional correlations have to be incorporated 
in order to properly describe the considered nuclei. The small energy dif- 
ferences observed between the coexisting minima indicate that many body 
effects beyond the mean field, like the restoration of the rotational symmetry 
and/or quadrupole fluctuations, may change the conclusions extracted from 
the raw Hartree-Fock-Bogoliubov approximation. The restoration of the ro- 
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tational symmetry leads to an energy gain (the so called rotational energy 
correction) which is proportional to the quadrupole deformation of the intrin- 
sic state and ranges from zero (spherical intrinsic state) to a few MeV for 
typical well deformed configurations in the region of the nuclear chart con- 
sidered in this study. Therefore, the consideration of this effect might play a 
key role for a more qualitative and quantitative description of the nuclei we 
are interested here. In addition, it is very important to consider the effect of 
quadrupole mixing [35-37,39,59,63]) in those cases where the angular momen- 
tum projected potential energy surfaces show important topological changes 
compared with the corresponding mean field surfaces. 

The theoretical background for the restoration of rotational symmetry is treated 
extensively in the literature (see, for example, Refs. [66-71]) and also very good 
books and papers concerning the treatment of correlations beyond mean field 
in the framework of the Generator Coordinate Method (GCM) are available 
(see, for example, Refs. [63,66,72]). Thus we are not going to dwell on these 
details here. In spite of this we believe that it is useful to present a short out- 
look of our procedure to carry out both angular momentum and configuration 
mixing in the case of a density dependent interaction like the Gogny force. In 
addition, the issue of which density dependence has to be used in these calcu- 
lations beyond mean field will be discussed in detail and strong arguments in 
favor of the prescription used will be given. 

2.1 Angular momentum projection with density dependent interactions. 

Let us consider a set of mean field states | (p(q)) depending on the parameters 
Q — {Q20, Q22, Q40, ■ ■ •} which define the corresponding mean field configurations. 
From this set we can build another set of states where the rotational symmetry 
is restored [66-71] 

\^iM{q))=Y.9 I K^P I M K \m)- (1) 

K 

The operator P[f K is the angular momentum projection operator [73] given 

by 

Pmk = I dnD* K {n)R(n) (2) 

with il representing the set of the three Euler angles (a,/3, 7), 2>^(f2) is the 
well known Wigner function [74] and R(Q,) = e ~ taJz e~' l/3Jy e~' l ^ Jz is the rotation 
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operator. The energy E T (q) is given by 





H 


*JAf(g)) 




ViM(q)) 



(3) 



The above procedure is based on the assumption that the hamiltonian is ro- 
tational invariant and therefore the intrinsic energy (tp(q)\H\(p(q)) is indepen- 
dent of the orientation of the intrinsic wave function \ip(q)), that is 

(<p($\H\<p($) = {<p@\R\n)HR(n)\<p@) (4) 



However, when dealing with density dependent (DD) forces the above assump- 
tion apparently breaks down as, in general, the density dependent term is not 
rotational invariant. This apparent paradox can be solved if density dependent 
hamiltonians are thought not as genuine hamiltonians but rather as devices 
to get an elaborated energy functional of the density. With this point of view 
the right hand side of Eq. (4) has to be treated as 

(^l^t^)^^^]^)!^) (5) 

where the interaction depends on the rotated density 

Ai,n(rO = (<pmR\to)ptfR(n)\<p($) (6) 



with the density operator defined in the usual way 

(r)=X>(r-r1) (7) 



Before proving that Eq. (4) holds for density dependent interactions let us 
present some basic results. First, we introduce the 3x3 unitary rotation 
matrix 71(H) = TZ z (a)TZ y (l3)TZ z (^) that appears when eigenstates | r) of the 
position operator are rotated 

R(n) | r) =| Kr). (8) 



This is the same matrix that appears when vector operators (like the position 
operator) are rotated 

R?(Q)?R(Q) — TZ(Q)r (9) 



Now we have to consider the role played by the density operator. Its mean 
value (or overlap) appears in the density dependent part of the hamiltonian 



6 



and therefore the "parameter" r of Eq. (7) has to be treated as the position 
operators in that case. As a consequence the action of the rotation operator 
on r has to be considered 

A A 

R\n)p(r)R(Q) =p(Jl(tt)r) = ^ S(TZ(Q)f - f<) = ^ 6(r - H\Q)ri)(10) 

i=i i=i 

as well as its action on the internal coordinates r*j 

A A 

R\Q)p(r)R(Q) = Y j 5{r-TZ{n)r l ) = ]T 5(R,\n)r - n) = p(7^(fi)f).(ll) 
i=i i=i 

With these results we have for the rotated density of Eq. (6) that 

ti(n)p n , n (r)R(n) = (^|Et(0)p(K(fi)r^(fl)|^-)) = (^|p(r)|^) = p{f) 

The first identity is the result of applying Eq. (10) whereas the second is 
the result of applying Eq. (11). The result thus obtained shows that even for 
density dependent forces Eq. (4) holds. 

In order to evaluate the projected energy we need to know the Hamilto- 
nian overlaps (ip \ Iv(Q) HR(Q') | p>). For density independent and rota- 
tional invariant interactions the hamiltonian overlap is simply given by (</? | 
HR(Q' — Q) | (p). For density dependent interactions this quantity is not 
well defined and therefore we have to look for a prescription for the density 
dependent term. The most popular prescription so far is to use the density 





R^{Q)p{r)R{Q') 






W(n)R(n>) 





(12) 



in the density dependent part of the interaction. Several arguments to show 
the credibility of this prescription have been given. For example, for those 
Skyrme forces, with a linear dependence in the density, that can be cast as a 
three body force the hamiltonian overlap can be computed without ambiguity 
and the resulting density dependent term depends on the density of Eq. (12). 
Arguments based on the extended Wick theorem (the fundamental tool to 
compute the overlap of operators between different product wave functions) 
are also very popular. Apart from the previous arguments, it is mandatory 
for the prescription to be used that it fulfills the following two requirements: 
First, it should not carry angular momentum or, in other words, the following 
identity should be fulfilled 

(V hM M | H | ^>i 2 mM) = S hh S Ml M 2 (^hMM | H | *i lMl (q))- ( 13 ) 
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Second, the prescription must lead to a real projected energy (notice that, in 
general, the density of Eq. (12) is a complex quantity). 

The first requirement (Eq. (13)) is a direct consequence of the property 



Paxv ( r ) = ; ^ , , tttt ( 14) 









Ho) 







Ho) 


p(TV{0)r)R{n' - 








r(w - n) 


¥>(?)> 



which is easily obtained using Eq. (11). Using this result we obtain 

R\n)p n ^,(r}R(n) = Po,n'-n(?)- (15) 

This property allows to write the numerator of the projected energy as 

/ dndQ'V^WV^in'X^q) | #(fi)jf [p^ n ,(r)} R(Q') \ V (q)) = 
J dmn'V^^V 1 ^^)^ | H [po^-atf] R(n> - fi) | <p(q)) = 
E / ^IWO^mUO / dQV^ K2 (Q)( V (q) \H[p n (r)]R({l) \ <p{q)) = 

171 J J 

^Y^.mAa/ dnV^m^q) \H[p Q (r)]R(n) I v^(g)) 

with 

Pn(r) = p ,n(0 

given by Eq. (12). 

The projected energy Eq. (3) can now be written as 

j Ekk' g£(g)&(g) J dnv^ K ,(n)h(n) 

[q) Ekk' gSWAti I dQV^ K ,(Q)n(Q) [ } 

with 

h(Sl) = (<p(q) | H [p n (r)} R(Sl) | <p@) (17) 

and 

n(Q) = (tp(q) | i?(fi) | cp(q)). (18) 
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This result shows that the energy is a quantity independent of the quantum 
number M and therefore of the orientation of the laboratory reference frame. 

In order to ensure the reality of the projected energy [67-69] we have to show 
that the property Q 

h*(Q) = h(-Q) (19) 



holds even for density dependent forces. Here we would like to remark again 
that the mixed density of Eq. (12) is not in general a real quantity and is 
not necessarily positive definite. The complex conjugate density in the density 
dependent term is given by 





R ] {n)p{r) 






i?t(0) 





_ (<p% I p(1l\tt)r)R(-tt) | y(g)) 
&(q) | R(-a) | <p($) 

(20) 

Using this result we have 



R j (m[pUr)]\ V (q)) 



R\n)H \p_ n (Ttf(tt)rj\ R(Q)R(-Q) \ ip(q)) 



(21) 



Now, using 

R\Q)p_ n (n\Q)f)R(n) = p_ n (f), 



(22) 



we finally obtain 



p-o\r, 



R(-n) | tp(q)) = h{-Q). 



(23) 



2 The same property holds for the norm overlap n(Q) and also for the neutron and 
proton overlap iV(fi) = (<p(q) \ NR{ty \ <f(q)) and = (ip{q) \ ZR{n) | ip{q)). 

Here and in the following we have mainly concentrated on the hamiltonian overlap 
h{p). 
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Besides the fact that with the prescription of Eq. (12) one recovers for the 
energy an expression similar to the one already known for density independent 
forces (Eq. (16)), it is also important to note that when the intrinsic wave 
function is strongly deformed, and the Kamlah expansion can be used to 
obtain an approximate expression for the projected energy (cranking model), 
a prescription like the one of Eq. (12) yields the correct expression for the 
angular velocity u including the "rearrangement" term [75,76]. 



2.2 Restriction to axially symmetric intrinsic wave functions 



For computational reasons, we restricted ourselves to axially symmetric (K = 
0) configurations. In this way two of the three integrals on the Euler angles 
can be performed analytically reducing by at least two orders of magnitude 
the computational burden. The axially symmetric Hartree-Fock-Bogoliubov 
(HFB) wave functions | <p(q2o)) were obtained in mean field calculations with 
the axial quadrupole moment Q20 — z 2 ~ 5 {% 2 + y 2 ) as constraining operator 
[fj Besides the axial symmetry, we have imposed as selfconsistent symmetries 
the parity, the e~ inJv symmetry and time-reversal. 

Due to the K = restriction in the HFB states | ^(920)) (i- e -; Jz | <£>(<?2o)) = 0) 
we obtain 



Pn(r) 



<¥>(?20) 


p(r)R(n) 




p(r) 


e -iaJz e - 


-i/3Jy 


¥>(?2o)) 






¥>(?2o)) (v(Qw) 


g T>(3J y 


V(?20)> 







g ijjj y 


V(?20)> 




g ifjj y 


¥>(?!»)) 



p p (Kl(a)r) (24) 



where TZ z (a) is the rotation matrix along the z axis. In the same way 



h(Q) = (<p( q20 ) I H [p n (fO] R(n) I ^(g 2 o)> 



(v(feo) 

M?2o) 



p /3 (^l(a)r) 



air 



-iaj z 



¥>(?2o)) 



r 



(25) 



Using the selfconsistent symmetry e l7rJi/ in the HFB wave functions | f{q2o)) 
and the identity e t/3Jy = e~ inJz e~ il3Jv e l7Tjz it can be easily shown that pg(r) and 

3 The definition used for the intrinsic quadrupole moment is a factor \ smaller 
than the usual definition for this quantity. 
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h{(3) are real quantities that also satisfy pg(r) = p n _Jr) and h(/3) = h{it — (3). 
Using these properties and — (3) = ( — ) 7 rfoo(/^) the integration interval on 
(3 can be reduced to [0, 7r/2] and it can be shown that the integrals are only 
different from zero when / is even (notice that the selfconsistent symmetry 
e -nrj y f orces positive parity intrinsic states). Therefore, the projected energy 
is a well defined quantity for even I and for odd / becomes an indeterminacy. 
The projected energy can finally be written, for each (^-configuration in the 
following form 



rt, v A/n /o 5 ^sin(/5)^(/3)(v9(g 20 ) \ H \p (?)} e'^y \ <p(q 20 )) 
E (q w )=A(I)- 



gdf3sm((3)d I M(( P (q 20 ) | e-Wy \ V (q 20 )) 

_ A{I) gdpsm(WMHP) (26) 

g dPsm(P)d^(P)n(P) 

where A(J) = \ (l + (-) J ) has been introduced to recall that for odd / the 
projected energy is an indeterminacy. In the above expression for the energy, 
the density to be used in the density dependent part of the interaction is given 
by 



M?2(>) 


p(f)e- il3J y 


y?teo)> 






¥>(?2o)) 



m = /. ' ?'„.,„ , " (27) 



Finally we would like to mention that in order to account for the fact that the 
mean value of the particle's number operator usually differs from the nucleus' 
proton and neutron numbers, we followed the usual recipe (see for example 
[69,72]) and replaced h{(3) by ti {(3) = h{(3)-X z {Z{(3) - Z )-\ N (N(f3) - N ) 
where Z{(3) = ( V (q 20 ) \ ZeT* 3 * \ <p(q 20 )}, N((3) = ( V (q 20 ) \ Ne'^ \ V (q 20 )} 
and \z and A^v are chemical potentials for protons and neutrons, respectively. 



2.3 Angular momentum projected configuration mixing with density depen- 
dent interactions. 



Once we have the angular momentum projected potential energy surfaces 
(AMPPES), the next step is to carry out configuration mixing (AMPGCM). 
In this case the ansatz for the wave function is 



*im{<t)) = I dqf'^q) | * IM (q)) (28) 
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where the wave functions | ^im{q)) are given by Eq. (1). The amplitudes 
are solutions of the Hill-Wheller (HW) equation [77,66]. As it was 
mentioned before, for the moment, we restricted ourselves to axially symmetric 
(K = 0) configurations and in this case the amplitudes / / ' <J (q , 2o) are found by 
solving the HW equation 

/ dq 20 (H I {q 20 , 4) " E I ^ I (q 20 , qj) f^{q 20 ) = (29) 
with the kernels H I (q2o, q'20) and N I {q2Q, q'20) defined, for even spin values, as 

2 

W z (g 20 ,<4) = ( 2/ + !) / ^sin(/3)4*(/3)(^( g20 ) | # [p? CM (r)] e"^ | ¥>(rio)>(30) 



and 

TT 

2 

A^(g 20 ,4) = (2/ + 1) / rf/3sin(/3)4*(/3)(vp(g 2 o) I e"^ | ^(g 20 )> (31) 



The previous results for H T (q 2 o, q' 20 ) an d A^ 7 (q 2 o, <7 20 ) can be found along the 
same lines described in the previous section. In Eq. (30) the density 
is given by 

-jfiCM 1 -a _ <^(g2o) I pe~ i/3Jw I <p(q' 20 )} 



which is the generalization of Eq. (27) for the density dependent part of the 
interaction in the framework of the configuration mixing calculation. As in 
the case of angular momentum projection we also replaced H (j)p CM (r)) by 

H (pp CM (r)) - \z(z - Zo) - \ N (N - A ) where \ z and are chemical 
potentials for protons and neutrons, respectively. 

The first step [66] in the solution of the HW equation Eq. (29) is to diagonalize 
the norm kernel 

J dq'^N 1 (q 20 , q^uliq'w) = n^(g 20 )- 

The norm eigenvalues n{ with zero values are subsequently removed (i.e, lin- 
early dependent states are removed from the basis) for a proper definition 
of the collective image of the kernel Ti. I (q 2 Q, 9 2 o) which is given by 

4 In practical computations eigenvalues smaller than a given threshold e should be 
removed to ensure the numerical stability of the solution of the HW equation. 
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H k£W = ri I dq20 I d V I C{q2Q)H I {q2o,q 2 Q)u I k >{q 2{) ) (33) 



ln ky n k > 

and the diagonalization in the collective subspace 3c 



EH^W/^E^gl' (34) 



I,a _ 

A± kk' 

kk' 

gives us the energy E I,cr for each spin value not only for the ground state 
(a = 1) but also for excited states (<r = 2,3,4, . . .) that, with the considered 
set of generating wave functions, could correspond to states with a different 
deformation from the one of the ground state and/or to quadrupole vibrational 
states. 

Using the eigenfunctions i^feo) of the norm kernel N I {q 2 Q ) q'20) an d the am- 
plitudes g]f (see Eq. (34)) we can compute [66] the amplitudes f I,<T (q2o) 

/ /lff teo) = E^teo) (35) 
k yn^. 

and the so called collective wave functions G I,a (q 2 o) 

G / ' CT (g 20 )=E^feo). (36) 

k 

These collective wave functions are orthonormal and therefore their module 
squared can be interpreted as a probability amplitude. 

Finally, from the knowledge of the amplitudes f I,a {q2o), we can compute the 
reduced B(E2) transition probabilities and the spectroscopic quadrupole mo- 
ments Q spec (I, a). This is one of the main motivations for carrying out a con- 
figuration mixing calculation of angular momentum projected wave functions 
in the case of Gogny [36,37,39,62,59] and Skyrme [63,35] forces, since both 
interactions allow the use of full configuration spaces and then one is able to 
compute transition probabilities and spectroscopic quadrupole moments with- 
out effective charges. In the framework of the AMPGCM the B(E2) transi- 
tion probability between the states (ij,<7j) and (If,af) is expressed as (see 
Appendix A for more details) 



B(E2,I m ^I f a f ) = ^ n (37) 
J dq2o : idq 20 jf I i> a f*(q 20 j)(Ifq2oj \ \ Q2 \ \ AftM,*) /^(ftio,;) 
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and the spectroscopic quadrupole moment for the state (I > 2, a) is given by 



16tt 12 1 

«*"<'• 'Wir (,„_,! (38) 

x J dq 20:i dq 20 jf I > a *(q20,i)(Iq20,i \ \ Q2 \ \ Iqw,f) / V (?20,/) 



2.4 Details on the calculation 



The intrinsic wave functions | <£>(<?2o)) were obtained as the solutions of the 
constrained axially symmetric HFB equations with the constraint on the mass 
quadrupole moment. The HFB creation and annihilation operators were ex- 
panded in an axially symmetric Harmonic Oscillator (HO) basis including ten 
major shells (i.e., Nsheii — 10). The two oscillator length parameters of the 
basis b± and b z were chosen to be always equal in order to keep the basis 
closed under rotations [78,79] (this is also the reason to include full HO major 
shells in the basis). The same oscillator length was used for all the quadrupole 
deformations considered in order to avoid completeness problems in the GCM 
calculations [78]. 

In the HFB calculations the two body kinetic energy correction was fully 
taken into account both in the energy and in the minimization procedure. 
This term is specially important in the nuclei considered due to their small 
mass number. Concerning the Coulomb interaction, both the exchange and 
pairing parts were not taken into account as they increase the computational 
burden by at least an order of magnitude. However, we computed the Coulomb 
exchange energy in the Slater approximation and added this quantity at the 
end of the calculations in a perturbative fashion. 

As self consistent symmetries we kept, in addition to the axial symmetry, the 
parity (no octupole mixing is allowed), the e~ tnJy symmetry and time reversal. 

For the computation of the matrix elements of the rotation operator in the 
HO basis we have used the results of [80]. 

For the evaluation of overlaps of one or two body operators between different 
HFB wave functions we have used the extended Wick's theorem [81]. In order 
to determine the sign of the norm overlaps we have followed the procedure 
proposed in Ref. [82]. 
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Fig. 1. On the left hand side, the ground state HFB energies for 40 Mg are plotted 



as functions of the HO length b Q (fej 



b ) for the Ngheii = 7, 8, 9, 10, 11, 12, 13 



and 17 bases. The curve corresponding to the basis used in the present work 
{Nsheii = 10) is plotted as a dashed line. On the right hand side, lower panel, 
the HFB energies computed with b Q = 2.1 fm and Nsheii = 10 and 11 are plotted 
as a function of the quadrupole moment. In the upper panel, the energy difference 
between the N 'shell = 10 and N shell = 11 calculations is plotted as a function of the 
quadrupole moment. 



3 DISCUSSION OF THE RESULTS. 



3.1 Mean field approximation for Magnesium isotopes. 



A few words concerning the convergence of our calculations with the size of the 
HO basis are in order here specially taking into account that we have to deal 
with the dripline nucleus 40 Mg. One should note that due to the proximity of 
the two neutron dripline, the full HFB approximation must be used [83,9] and 
absolute convergence for the binding energy can only be found for HO basis 
with a very large number of shells. At the mean field level such a hard task is 
still feasible. Even in the case we were interested in a single g2o configuration, 
angular momentum projected calculations with very large Nsheii can also be 
performed. The main reason why we can not consider so large Nsheii m the 
present study is obvious: the enormous amount (typically of the order of one 
thousand) of angular momentum projected hamiltonian overlaps to be com- 
puted in the configuration mixing calculation make the problem intractable. 
In 40 Mg (using the 520 range —2.06 < ^20 < 3.06, the mesh Ag 2 o — 0.16 and 
16 points in the angular momentum projection) if the basis is increased from 
N shell = 10 to N shell = 11 the computational time required to evaluate all 
the projected kernels increases by a factor of 10 (from 3 days to one month 
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in a typical workstation). Here it is worth to remark that besides the tech- 
nical difficulties in the configuration mixing calculation, our force is a finite 
range force and its mathematical structure leads to very complicated matrix 
elements whose evaluation is very time consuming. 

On the other hand, as the absolute value of the binding energy does not affect 
very much the collective motion (it is only affected by the shape of the energy 
landscape) we can select a reasonable value of N S h e u for which the energy 
landscape is well converged (i.e. its shape remains almost unchanged in the 
regions of physical significance when the value of N S heii is increased). As we 
will see in the next paragraph Nsheii = 10 is more than enough for the present 
study. 

On the left hand side of Fig. 1 we show the ground state HFB energy for 40 Mg 
as a function of the oscillator length b Q (b± = b z = b Q ) for different values of 
N shell (the Coulomb exchange energy is not included). As expected the curves 
become more and more flat for increasing values of Nsheii ■ Already the Nsheii = 
17 basis can be considered as a reasonable approximation for an infinite basis 
in this nucleus as the dependence of the energy on the oscillator length is very 
weak for a wide range of b Q values. For N She u = 10 a minimum in the energy 
curve is obtained for b Q = 2.1. Using the minima of the corresponding energy 
curves for N Sh eii = 10 and N S hdi = 17 we get E NshM=w - E NshM=17 = 568.37 
keV. In 38 Mg, the same analysis has been carried out and the overestimation 
of the energy is 543.29 keV. As a consequence, for the two neutron separation 
energy S 2N it is obtained that | S 2N ,N ShM =io{ m Mg) - S 2N ,N ShM =n{ m Mg) | = 
25.12 keV. Also a very good agreement is found, in 40 Mg, between our values 
for the proton and neutron f3 2 deformation parameters [fi 2 , f3 2 ) = (0.38, 0.31) 
and the ones (0.35, 0.28) of the calculation in coordinate space of Ref. [29]. 
This can be understood from the fact that our basis roughly corresponds to 
R box ~ 9.4/m while it was found in Ref. [29] that the deformation parameters 
remain practically unchanged for R box > 7fm in 40 Mg. 

On the right hand side of Fig. 1 we show, in the lower panel, the energy 
landscapes of 40 Mg as a function of the quadrupole moment for the calculations 
with N shell = 10 and N S heii = H- In the upper panel we have represented 
the energy difference between both calculations. From this plots we observe 
how in the region between -1.5 b and 1.5 b the energy landscape does not 
change much when the basis is increased. As we will see in another section 
this range of quadrupole deformation is the one where the collective dynamics 
is concentrated and therefore it is not expected to find significant differences 
between the calculations with N S heii = 10 and N S heii = 11- 

Since the main interest of the present study is focussed on physical quantities 
like S 2 n values, rotational energy corrections, excitation energies, etc., and 
these quantities do not change very much with the number of HO shells, we 
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Fig. 2. Mean field potential energy surfaces for Magnesium isotopes as functions 
of the axially symmetric quadrupole moment. These curves have been shifted to 
accommodate them in a single plot. The corresponding energy shifts in panel a) are 
-85, -65, -50, -35, -25 MeV for 2 °- 28 Mg, while in panel b) they are -10, -7, -7, -8, 
-15, -20 MeV for 3 °- 40 Mg. 

conclude that our mean field results with Nsheii = 10 can be regarded as 
a reasonable choice. Even more, the zero point rotational energy correction 
Erot = Ehfb — E I=0 is pa 3 MeV for the ground state configuration in 40 Mg, 
i.e, the zero point rotational energy correction is 5.3 times larger than the 
relative energy difference E NshM=w — E NshM=17 , and also remains practically 
unchanged for increasing values of N S heii- 



In Fig. 2 we show the mean field potential energy surfaces (MFPES) as a 
function of the axially symmetric quadrupole moment q 2 o for the even-even 
Magnesium isotopes 20 ~ 40 Mg. The MFPES shown do not include the Coulomb 
exchange energy and they have been shifted to accommodate them in a single 
plot (see caption). As one can see, with the exception of 20 Mg, the MFPES are 
very flat around the corresponding minima indicating that further correlations 
could, in some cases, change the conclusions obtained at the mean field level. 

The nucleus 20 Mg shows a well denned spherical minimum which is a conse- 
quence of the N=8 neutron shell closure. On the other hand, both 22 Mg and 
24 Mg are prolate deformed in their ground states. In 22 Mg the ground state 
corresponds to g 2 o — 0.55(/3 2 = 0.40) and an oblate local minimum also ap- 
pears at g 20 = —0.26 (f3 2 — —0.17) with an excitation energy of 1.65 MeV. 
In the case of 24 Mg the ground state corresponds to g 20 = 0.66 (/3 2 = 0.43) 
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and another local minimum is found at q 2 o — —0.36 (fa — —0.22) with an 
excitation energy of 3.86 MeV. The only oblate isotope in this chain is the 
nucleus 26 Mg whose ground state is located at q 20 — —0.46 (fa = —0.26). A 
prolate isomeric state is also found at q 20 = 0.56 (fa = 0.32) with an excitation 
energy of 707 keV with respect to the oblate ground state. On the other hand, 
the nuclei 28 Mg, 30 Mg and 32 Mg show spherical ground states. The MFPES 
of both 28 ' 30 Mg are particularly flat around their spherical ground states. In 
the nucleus 32 Mg we obtain a prolate shoulder at q 20 — 0.86 (fa = 0.36) in 
the MFPES with 1.86 MeV of excitation energy with respect to the spherical 
ground state. 

From 34 Mg to the dripline isotope 40 Mg, prolate deformed ground states are 
found. The ground state deformations are q 20 — 0.96 (fa = 0.36), q 20 = 0.96 
(fa = 0.33), q 20 = 1.06 (fa = 0.33) and q 20 = 1.16 (fa = 0.33) in 34 Mg, 
36 Mg, 38 Mg and 40 Mg, respectively. Besides the fact that prolate configurations 
become dominant in all these nuclei, one should note that close-lying oblate 
isomeric states are found in the MFPES of 36 Mg, 38 Mg and 40 Mg at q 2 o = —0.46 
(fa = -0.15), q 20 = -0.66 (fa = -0.20) and q 20 = -0.86 (fa = -0.25) with 
excitation energies of 2.98 MeV, 2.51 MeV and 1.38 MeV with respect to 
prolate ground states. 

In Fig. 3 the proton and neutron particle-particle correlation energies ( defined 
as — E pp = ^Tr (A/c*)) are plotted as a function of the quadrupole deformation 
for all the isotopes considered. The evolution of the particle-particle correlation 
energies is well correlated with the structures found in the MFPES. Non-zero 
proton pairing correlations are found in all the spherical or oblate minima. 
In addition, sizeable neutron pairing correlations are found in 22 ' 24 Mg and in 
36 ' 38 ' 40 Mg for the spherical and oblate minima. Vanishing proton pairing cor- 
relations are found in the prolate side starting at q 2 o = 0.5b in all the isotopes. 
The range of quadrupole moments for which the proton pairing correlations 
vanish increases with the neutron number. An expected result is that neutron 
pairing correlations vanish in the mean field spherical ground states of both 
20 Mg and 32 Mg as a consequence of the N=8 and N=20 shell closures. However, 
tangible neutron pairing correlations are found at the spherical configuration 
of 40 Mg pointing to the erosion of the N=28 spherical shell closure already 
at the mean field level. On the other hand, the lowering of the level density 
around the ground state of this nucleus (q 2 Q ~ 16) leads to vanishing proton 
and neutron pairing correlations. This agrees with the results of [60,61,38] and 
[59] which suggest that while the N=28 spherical shell closure disappears for 
some neutron rich nuclei a new deformed shell closure emerges on them. 

It should be stressed here that the unphysical collapse of pairing correlations, 
which is clearly visible from the results of Fig. 3, is one of the main drawbacks 
of our calculations at the present stage. One should also consider the dynamical 
pairing correlations in these nuclei and their coupling with the quadrupole 
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Fig. 3. In each panel, proton (thick dashed lines) and neutron (thick full lines) par- 
ticle-particle correlation energies — E w are depicted as a function of the quadrupole 
moment. Also the HFB energy curves are plotted as thin full lines. 
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degree of freedom in the scope of a theory beyond the mean field in order to 
treat in an equal footing both short and long range correlations. 

Using the absolute minima of the MFPES we have computed the two neutron 
separation energies S 2N = E MF (Z,N — 2) — E MF (Z,N) and the results are 
compared in Fig. 4 with the ones previously obtained in the framework of 
the mean field approximation with the Skyrme force SIII [29] and SLy4 [30] 
and with the available experimental values [84]. Our results agree quite well 
with those of SIII in all the nuclei studied. In the case of SLy4 the agreement 
is also rather good except for 22 ' 24 Mg and 36 Mg where the SLy4 results are 
bigger than ours. The three interactions predict that the nucleus 40 Mg is the 
last bound isotope of the chain. Also the three interactions predict a dip in the 
S2N of 34 Mg and therefore fail to reproduce the peak observed experimentally. 
The observed peak is usually explained as a consequence of a deformed ground 
state in 32 Mg which has a lower energy than the spherical configuration. As 
we will see later, the wrong behavior of the S 2 n of 34 Mg at the mean field 
level is cured when the extra correlation energy (coming from the interplay of 
the quadrupole fluctuations and the zero point rotational energy corrections) 
leading to a deformed ground state in 32 Mg is considered. Finally, another 
significant failure of both the Gogny and Skyrme SIII interactions is that of 
the predicted values of the two neutron separation energies of 22,24 Mg. As 
we will see later, considering the zero point rotational energy correction will 
improve that situation. 

In Fig. 5 we show the single particle energies for protons and neutrons in 
32 Mg. As in our mean field calculations we solve the full HFB equations, the 
only quantities that can be properly defined are the quasiparticle energies. 
However, in order to have the usual Nilsson like diagrams we have chosen to 
plot the eigenvalues of the mean field Hartree-Fock hamiltonian h = t + T as 
a function of the axially symmetric quadrupole deformation. Looking at the 
neutron single particle energies we observe how at q 2 o = 0.86 a couple of f 7 / 2 
orbitals cross the Fermi surface and become occupied. The occupancy of those 
orbitals leads to the appearance of the shoulder seen in the MFPES of 32 Mg. 
It is also interesting to notice that the almost full occupancy of the d 5 / 2 proton 
orbital favors oblate shapes. In the same way, the full occupancy of the d 5 / 2 
neutron orbital in 26 Mg favors an oblate ground state for this nucleus. 



3.2 Correlations beyond the mean field : Angular momentum projection and 
configuration mixing. 

The main outcome of our angular momentum projected calculations is pre- 
sented in Fig. 6, where the angular momentum projected potential energy sur- 
faces (AMPPES) are shown for the nuclei 20_40 Mg. The corresponding mean 
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Fig. 4. The two neutron separation energies obtained in the present calculations are 
compared with the experimental data taken from Ref. [84] and with the mean field 
results obtained with the Skyrme interactions SIII [29] and SLy4 [30]. 

field energy landscapes (dashed curves) are also included for comparison. 

The first noticeable fact that one can see is that, with the exception of the 
I w = + curves, several points around the spherical configuration have been 
omitted in the AMPPES for the 1^0 curves. The reason for such omission can 
be understood by looking at Fig. 7 where the behavior of the projected norm 
A/" / (g 2 o, 920) = (0(920) I Pqo I 0(92o)) as a function of q 2 o for the nucleus 32 Mg 
is shown. From this figure it becomes clear that the omitted points, in this 
and the other nuclei, correspond to intrinsic configurations with a very small 
value of the projected norm A/" 7 (920, 920) • In those situations the projected 
energy is the quotient of two very small quantities and therefore its evaluation 
is affected by numerical inaccuracies that lead to erratic values (deviations 
from the smooth trend can be as large as a few MeV) . Fortunately, due to the 
smallness of the projected norm, these points can be safely omitted since they 
do not play any role in the configuration mixing calculations [37,39,62,63,35,59] 
to be discussed later on. We also observe that no energy gain is obtained for 
the spherical configurations and I n = + . From a physical point of view this 
is the expected behavior since these spherical configurations are already pure 
+ states as can also be seeing in Fig. 7 where we have A/" /=0 (0, 0) = I. 

Comparing the MFPES and the I n = 0+ AMPPES of the nuclei 20 Mg, 28 Mg, 
30 Mg and 32 Mg one can see that while at the mean field level spherical ground 
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Fig. 5. Proton and neutron single particle diagrams as a function of the quadrupole 
deformation for the nucleus 32 Mg. Full (dashed) lines correspond to positive (neg- 
ative) parity levels. The Fermi levels are represented as thick dotted lines. 



states are obtained, once we carry out angular momentum projection two 
minima, one prolate and the other oblate, appear. These minima are very 
close in energy and the difference amounts to around 40, 80, and 49 keV 
in 20 Mg, 28 Mg, 30 Mg, respectively, and 637 keV in 32 Mg. These minima are 
separated by a spherical barrier which is around 1, 3, and 2.6 MeV in 20 Mg, 
28 Mg, 30 Mg, respectively, and 1.6 MeV in 32 Mg. 

It is noticeable to observe how the ground state of 32 Mg becomes deformed 
after the inclusion of the rotational energy correction. The intrinsic configu- 
ration corresponding to the shoulder seen in the MFPES at q2o = 0.86 has a 
big correlation energy coming from the restoration of the angular momentum 
that is big enough as to overcome the energy difference with the spherical 
configuration. The deformed intrinsic configuration at ^20 — 0.86 corresponds 
to (see Fig. 5) a configuration in which two neutrons from the fan shell have 
crossed the Fermi surface. In the shell model language this is just a config- 
uration where a pair of neutrons have been promoted from the sd to the pf 
shell and this is the physical mechanism invoked [41,43] by the shell model 
practitioners to explain deformation in this nucleus. 

All the other Magnesium isotopes, with the exception made of 26 Mg are prolate 
deformed at I* = + . In fact, for increasing spin values the prolate minima 
become deeper than the oblate ones or the oblate minima are washed out. 
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Fig. 6. Angular momentum projected potential energy surfaces (full lines) for the 
nuclei 20_40 Mg and for the spin values I n = + , 2 + , 4 + , 6 + as functions of the axially 
symmetric quadrupole moment q2o ■ The mean field potential energy surfaces are also 
plotted as dashed lines. In the nucleus 40 Mg we have also included (dotted lines) 
the projected results corresponding to the calculation with N s h e u = 11. 

The nucleus 26 Mg, with its oblate intrinsic state at I n = + is the exception 
in all the Mg isotopes studied. As it was mentioned before, the responsible for 
such oblate minimum is the full occupancy of the neutron d 5 / 2 orbital which 
favors oblate deformations. The intrinsic state of the lowest lying I w = 2 + 
state remains oblate deformed but already at I n = 4 + the underlying intrinsic 
state becomes prolate deformed. The spectroscopic quadrupole moment of 
the lowest 2 + state in 26 Mg is known experimentally [26] to be -13.5 (20) fm 2 
indicating that this is a prolate deformed state. However, the experimental 
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Fig. 7. Projected norm J\f I (q2o,Q2o) as a function of q 2 o for the nucleus 32 Mg. For 
more details see the main text. 

low excitation energy of the state in this nucleus (3.588 MeV, versus the 
6.432 MeV for the same quantity in 24 Mg) is an indication of a strong shape 
coexistence between the prolate and oblate solutions. 

From the above discussed results we realize that the AMPPES for some nuclei 
and some spin values show the phenomenon of shape coexistence and therefore, 
a configuration mixing calculation is needed to disentangle the structure of 
those states. Even in those situations where the AMPPES show a relatively 
well pronounced minimum it is always interesting to check its stability in the 
framework of the configuration mixing calculation. The reason is that not 
only the AMPPES but also the collective inertia has to be considered in the 
framework of a dynamical treatment in order to determine the stability of a 
given solution. With this fact in mind, we have carried out Angular Momentum 
Projected Generator Coordinate Method (AMPGCM) calculations along the 
lines described in subsection 2.3 using the intrinsic axial quadrupole moment 
q 2 o as generating coordinate. These configuration mixing calculations have 
been performed with a mesh Aq 2 o = 10/m 2 which was tested to be accurate 
enough to describe, at least, the low-lying spectrum we are interested in this 
study. 

In Figs. 8, 9, and 10 we show the ground state (a = 1) collective wave func- 
tions squared | G ,/ '°" =1 (g2o) f fo r all the Magnesium isotopes studied in this 
paper up to I w = 6 + . We also plotted in each panel the AMPPES for the 
corresponding spin value. It is worth pointing out that from the position of 
the tails of the wave functions relative to the projected energies (see figure 
caption) we can read the energy gain due to the quadrupole fluctuations. In 
order to understand in a more quantitative way these collective wave functions 
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Fig. 8. Collective wave functions squared for the ground states (a = 1) and the spin 
values I 77 = + ,2 + ,4 + ,6 + of the nuclei 20 ' 22 ' 24 Mg. The corresponding projected 
energy curve is also plotted for each spin value. The y — axis scales are in energy 
units and always span an energy interval of 15 MeV (minor ticks are 0.5 MeV apart). 
The collective wave functions have also been plotted against the energy scale after 
proper scaling and shifting, that is, the quantity E I,a + 15 x | G I,a=1 {q2o) \ is the 
one actually plotted. With this choice of scales we can read from the figure the 
energy gain due to the quadrupole fluctuations by considering the position of the 
wave functions' tail relative to the projected curve. 

it is convenient to analyze [37,39,62,59] the averages 

120= J dffcol G^Ofco) | 2 g 20 (39) 

that give a measure of the deformation of the underlying intrinsic states. 

The Of wave function for 20 Mg shows a great admixture between the prolate 
and the oblate minima found in the AMPPES. In fact, prolate and oblate 
configurations have practically the same weight and therefore the ground state 
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Fig. 9. The same as Fig. 8 but for the nuclei i6 ' i8 '' 30 Mg. The y axis scales span in 
this case an energy interval ofl3 MeV. 

of this nucleus is spherical on the average. The energy gain due to quadrupole 
fluctuations is only 582 keV. As a result, the deformation effects previously 
found in the AMPPES are not stable once quadrupole fluctuations are taken 
into account and N=8 remains, on the average, as a spherical magic number. 
On the other hand for spin values I n > 2 + the ground state collective wave 
functions are almost inside the prolate wells, with the average deformations 
q 2 o = 0.376, 0.536 and 0.646 for the states 2^,4^, and 6f respectively. 

For 22 Mg and 24 Mg the ground state collective wave functions are well inside 
the prolate wells. This clearly shows on the one hand, the stability of the 
deformation effects found in the corresponding AMPPES and on the other 
that both systems are dominated by prolate deformations in the considered 
spin range. The average deformations are 0.416, 0.586, 0.626 and 0.636 for the 
states Of, 2f, Af, and 6f in 22 Mg, while the corresponding values in 24 Mg are 
0.596, 0.696, 0.706, and 0.766. 

In the nucleus 26 Mg, both the 0^ and the 2f states are slightly oblate deformed 
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Fig. 10. The same as Fig. 9 but for the nuclei 3 40 Mg. 

with g2o~ 0cr_1 = —0.176 and qi() 2,a ~ 1 = —0.166 while for the state 4f the col- 
lective wave function becomes prolate deformed (i.e. a band crossing takes 
place) with q 2 ^ a=1 = 0.396 and q 2 Q 6 ' a=1 = 0.726. In both 28 Mg and 30 Mg the 
ground state shows considerable mixing between the oblate and prolate con- 
figurations. Here for the spin values I n > 2 + , the collective wave functions are 
almost inside the prolate wells. The average deformations are 0.506, 0.606 and 
0.676 for 28 Mg for the states 2+, 4^, and Qf while for 30 Mg the corresponding 
deformations for the same spin values are 0.586, 0.966 and 1.096. 

The oblate character of the + and 2 + states in 26 Mg already obtained without 
configuration mixing is preserved when it is included (although in the latter 
case almost spherical configurations are obtained), in contradiction with the 
experimental result (prolate character) extracted from the value of the spec- 
troscopic quadrupole moment of the 2 + state (-13.5 (20) fm 2 ). However, our 
results predict a strong shape coexistence for those states in 26 Mg as well as 
for the + states of 28_30 Mg. A characteristic fingerprint of shape coexistence 
comes from the position of the 0^ excited state: it is expected to lie at a rather 
low excitation energy in those situations. Experimentally [25,26], the excita- 
tion energy of the 0^ state is known in 24 Mg (6.432 MeV), in 26 Mg (3.588 
MeV) and in 28 Mg (3.862 MeV). The sudden drop in the 0^" excitation energy 
in going from 24 Mg to 26 Mg is a clear indication of shape coexistence in the 
latter (and also in 28 Mg) nucleus. In our calculations, apart from the shape 
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of the collective wave functions, we get for those 0^ excitation energies the 
values 5.675, 2.592 and 3.700 MeV for 24 ~ 28 Mg, respectively, that £1X6 db clear 
manifestation of shape coexistence in the last two isotopes. The inclusion of 
dynamical pairing correlations may improve our description of 26 Mg. Dynam- 
ical pairing correlations for neutrons will not change much the oblate side due 
to the large single particle gap between the and the S1/2 orbitals (see Fig. 
5). However, in the prolate side they can bring into the wave function the 
fj/2 orbital which has a big correlation energy coming from the restoration 
of the rotational symmetry. As a consequence, the prolate side will gain more 
energy than the oblate one turning the ground state of 26 Mg from slightly 
oblate to prolate. Although the inclusion of dynamical pairing in the present 
calculations is rather cumbersome, work along this line is in progress. 

Now in 32 Mg, the Of collective wave function shows a significant mixing of the 
oblate and prolate configurations and as a consequence the deformation in the 
ground state is reduced from the 1.06 obtained taking the absolute minimum 
in the I n = + AMPPES to <?2o~° ,<T_1 = 0.436 once quadrupole fluctuations are 
taken into account via configuration mixing. This shows that the presence of a 
deformed ground state in this nucleus is the result of a subtle balance between 
the zero point corrections associated with the restoration of the rotational 
symmetry and the fluctuations in the collective parameters (in our case the 
axially symmetric quadrupole moment). On the other hand, the presence of 
a deformed ground state indicates that for this nucleus N=20 is no longer a 
magic number. The 2f, Af, and Qf wave functions in this nucleus are inside 
the prolate wells and the average deformations are 0.886, 1.016 and 1.086. 

From 34 Mg to 40 Mg, the ground state collective wave functions become strongly 
prolate. The dynamical deformations of the 0^ states are 0.796, 0.776, 0.796, 
and 1.176 respectively. In all these nuclei, the a — 1 collective wave functions 
are well inside the prolate wells in the whole spin range considered. The dy- 
namical deformation values for the states 2f, 4f and 6f are 1.06, 1.136, 1.196 
in 34 Mg, 0.996, 1.036, 1.056 in 36 Mg, 1.036, 1.106, 1.146 in 38 Mg and 1.236, 
1.256, and 1.266 in 40 Mg. The results show the stability of the deformation 
effects in Magnesium isotopes as we move towards the dripline. Moreover, the 
presence of a deformed ground state in 40 Mg also points towards the erosion of 
the N=28 shell closure. Contrary to the N=20 case, the erosion of the N=28 
shell closure already appears at the mean field level and therefore we can say 
that it is "stronger" than for the N=20 case. To summarize, it can be con- 
cluded that, while N=8 remains a spherical magic number in the Magnesium 
isotopic chain, both N=20 and N=28 spherical shell closures do not remain. 

Here we will make a few comments on the results obtained in the framework 
of the configuration mixing approach for the nucleus 40 Mg when the basis is 
increased from Nsheii = 10 to Nsmu = 11- The effect on the projected energy 
landscapes can be seen in the corresponding panel of Fig. 6 where the Nsheii = 
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Table 1 

Ground band spectroscopic quadrupole moments Q spec (I, a = 1) in efm 2 for = 
2+, 4+, 6+ in the nuclei 20 ~ 40 Mg. The Shell Model predictions from [44,46] are shown 
in brackets. 

11 curves are plotted as dotted lines. Although there are differences (specially 
at large absolute values of g 20 ) between the N She u = 10 and N S heii = 11 curves, 
these differences are almost independent of the considered spin. Therefore, 
we do not expect big changes in the excitation energies of the ground state 
rotational band as is the case: these excitation energies are, for all spin values, 
around 10 keV higher for N SheU = 11 than for N S heii = 10. As a consequence, 
the transition gamma ray energies remain unaltered by the increase of the 
basis size. On the other hand, the excitation energies (with respect to the 
true ground state) of the members of the excited rotational band decrease on 
the average by 50 keV and therefore, as in the previous case, the intraband 
gamma ray energies remain the same. From this results we can conclude that 
our calculations are well converged in terms of the basis size. 

The previous results suggest that all the considered nuclei are dominated by 
prolate deformations. This is clear from the results we have obtained for the 
average deformation q^ -1 and also from the negative values of the ground 
band spectroscopic quadrupole moments presented in Table 3.2. With the ex- 
ception of 26 Mg, our results for the spectroscopic quadrupole moments show 
a very good agreement with the Shell Model predictions of Ref . [44,46] (these 
predictions are shown in brackets) for the spectroscopic quadrupole moments 
of the 2f states. Experimentally, the 2 + spectroscopic quadrupole moment of 
24 Mg is -16.6 (6) e 2 fm 2 [26] and the one of 26 Mg is -13.5 (20) e 2 fm 2 [25]. 
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For 24 Mg we obtain a reasonable agreement with experiment (a 15 per cent 
discrepancy) whereas for 26 Mg the disagreement is strong. In the latter case, 
we have already traced back the disagreement to the strong shape coexis- 
tence obtained in our calculations and also to the effect of dynamical pairing 
correlations (not included in the present work). 

The values obtained for the quadrupole moments g^cT °f the intrinsic states 
can be used to classify in terms of bands each of the physical states provided 
by the AMPGCM approach. In Fig. 11 we have plotted the energies of the 
AMPGCM states E 1 ^ in a diagram of energy versus quadrupole moment. 
Each energy E I)G is placed at a q 2 § value corresponding to q~20 ■ In addition, 
we have plotted the AMPPES for / = to guide the eye. Although, in the case 
of the a — 1, many of the features observed in this figure have already been 
discussed, there we also show the results corresponding to the first excited 
states (a = 2) provided by the AMPGCM approach. Note that one of the 
main advantages of such representation is that the band structure of each 
nucleus can be observed at a glance. 

In Fig. 12 we compare the results for the AMPGCM two neutron separation 
energies S 2 n = E Q +(N — 2) — E +(N) with the corresponding mean field 
results (see subsection 3.1) and also with the available experimental values 
[84]. The AMPGCM binding energy is the sum of the mean field binding 
energy of the intrinsic state plus the energy gain due to the restoration of the 
rotational symmetry plus the energy gain due to the configuration mixing. 
Therefore, the differences in the two neutron separation energies obtained in 
the AMPGCM and the mean field are due to the later two contributions. An 
analysis of those contributions show that the rotational energy correction is 
the main responsible for the differences observed in the two neutron separation 
energies. The AMPGCM S 2 n energies differ substantially from the mean field 
ones in the nuclei 22 Mg and 34 Mg and are much closer to the experiment. On 
the other hand, it is worth to remark that the nucleus 40 Mg remains the last 
bound isotope in the chain in both theoretical approaches. 

In Fig. 13 the excitation energies of the 2f and 4+ states and the B(E2, Of —> 
2f) transition probabilities obtained in the framework of the AMPGCM are 
compared with the available experimental values [25,15,16,22,19,20,23] and 
also with the predictions of the Quantum Monte Carlo Shell Model [45] and the 
Shell Model [44,46]. Concerning the B(E2) transition probabilities we clearly 
see, from panels a) and b), that the agreement with the available experimental 
data is rather satisfactory and in most cases (with the exception of 30 Mg 
where our prediction appears a little bit too low) our results stay within the 
experimental error bars. On the other hand, our results are also consistent 
with the predictions of the Quantum Monte Carlo Shell Model [45] and the 
Shell Model [44,46]. Our results, although not as good as the SM or QMCSM 
ones, are very satisfactory considering that the parameters of the force have 
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Fig. 11. The AMPGCM energies E 1 * for a = 1 and a = 2 and F = 0+,2 + ,4+,6+ 
are plotted in an energy versus quadrupole moment diagram for the nuclei 2 °- 40 Mg. 
The quadrupole moment of each AMPGCM state is given by the average quadrupole 
moment q^Q ■ The AMPPES for I w = + are also plotted to guide the eye. 



not been fitted to the region and/or the physics of quadrupole collectivity 
and also that no effective charges have been used in our calculations of the 
transition probabilities. 

The calculated excitation energies for the 2f and 4^" states are plotted in pan- 
els c) and d). They reproduce quite well the experimental isotopic trend and 
are again consistent with the theoretical trends predicted by both the Shell 
Model [44,46] and the Quantum Monte Carlo Shell Model [45] . In most cases, 
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Fig. 12. The two neutron separation energies S^jv for Magnesium isotopes as ob- 
tained in the framework of the AMPGCM are compared with the corresponding 
mean field values and also with the available experimental values taken from [84]. 
The AMPGCM S 2N are defined as E Q + {N - 2) - E Q + (N). 

however, our values come too high as compared with the experiment. In our 
previous works [37,59] we also noticed the same behavior in some N 20 and 
N ?s 28 nuclei (our predictions were too high as compared with the experi- 
ment). Probably a proper treatment of some missing correlations will give the 
quenching factors we need for a much better agreement with the experiment. 
Although it is very difficult to assert before hand what are the missing correla- 
tions, the mixing of different K values with a full triaxial angular momentum 
projection and a beyond mean field treatment of the dynamical pairing fluc- 
tuations can be important ingredients for a more realistic description of the 
nuclei studied in this paper. Another source of discrepancy could be related 
to the fact that ours is a calculation of the Projection After Variation (PAV) 
type instead of the more complete Projection Before the Variation (PBV). 
Having all this in mind and the free parameter character of our calculation we 
conclude that our results for the excitation energies of the 2f and 4^ states 
show a rather satisfactory agreement with the available experimental and are 
consistent with other theoretical predictions. 



4 Conclusions. 

With the aim to describe the phenomenology of quadrupole deformation in 
light nuclei we have performed calculations with the Gogny force and in the 
framework of beyond mean field theories for the nuclei 20 ~ 40 Mg. First of all, the 
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Fig. 13. The excitation energies of the states 2f and 4+ provided by the AMPGCM 
and the B(E2,0f — > 2f ) transition probabilities in 20_40 Mg are compared with the 
available experimental data [25,15,16,22,19,20,23] and with the theoretical predic- 
tions of the Quantum Monte Carlo Shell Model [45] and the Shell Model [44,46]. 



results show the fundamental role played by the angular momentum projection 
for a proper description of the physics under study. The effect of projection 
in all the physical observables is so big that it can not be overlooked as it has 
been common practice in many previous calculations. In addition, we also find 
that the effect of configuration mixing is in most cases rather relevant. From 
our results we conclude that for the Magnesium isotopes the N=8 shell closure 
is preserved whereas for N=20 and N=28 deformed ground states appear in 
the calculations. The three isotopes from 36M9 to the drip line nucleus 40 Mg 
are predicted to be prolate deformed in their ground states. Concerning the 
excitation energies and B(E2) transition probabilities of the low- lying excited 
states we obtain a reasonable agreement with experiment. The agreement is 
not as good as the one obtained with other approaches like the Shell Model 
or the Quantum Monte Carlo Shell model. The reason for that probably lies 
in the fact that our treatment of the problem, although it contains the most 
important ingredients, is not as refined and complete as the one of the SM 
and/or the QMCSM. In addition, the interaction used has not been fitted 
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to this specific region of the periodic table. The later could be consider as a 
drawback of our calculations but we think it is the other way around and in 
fact is a manifestation of the strong predicting power of the Gogny interaction. 
It is rather satisfactory to obtain the results presented in this paper with 
an interaction that is also able to reproduce, for instance, the fission barrier 
heights of 240 Pu. 

Finally, the consistency of the prescription used for the density dependent 
term of the interaction in the present calculations beyond mean field has been 
discussed in great detail. 
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5 Appendix A: Calculation of transition probabilities. 

In this appendix we present the basic formulas for the computation of an- 
gular momentum projected transition probabilities in the framework of the 
AMPGCM. The starting point is the transformation property of the multi- 
pole operators under rotations 



Using the well known result for the product of two Wigner functions [74] as 
well as the definition of Eq. (2) for the angular momentum projection operator 
and the property 
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(41) 



we obtain after some algebra the result 



P I K , ' f M f Qx»PM l K t = (hMM | I f M f ) 



(42) 
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With the definition of the projected wave functions of Eq. (1) and the previous 
result we obtain 

(*i t M,m I I *W$> = {m t^ 1 IfMf) Mt II & II W ( 43 ) 



yfilf + 1 



with 



M-W Q*\\M) = {2Ii + 1 E If + 1) U It - x E R%£>7)<4(£) 
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KiKjVjx 
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7, A I f 



\u fi -K f 



cMV^m^qf) | Q v i2(fl) | ^)> (44) 



Using now the expression (28) we get 



x | dqdqff If ' af *(qf)(I f qf II <?a || /ifc)/^*). 



(45) 



Finally, the expression for the B(E\, l^i — > //07) transition probability is 
written as 



2/i + l 



IT 4- 

J dqdqff^mWf 11 qa 11 mf^m 



(46) 



In the present work we are interested in the calculation of transition probabil- 
ities for axially symmetric HFB states labelled by the quadrupole deformation 
q2o- Taking advantage of the axial symmetry of the intrinsic wave function as 
well as the selfconsistent symmetry e~ lnJy we can simplify the above expres- 
sions as follows. First we have 

(¥>(«»,/) I QxA( Q ) I ¥>(<&o,i)> = e^V(?2o,/) I Qv e "^* I V(«»,i)> (47) 



that leads to 
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20,i)) = 47T 5 Q _^5 Ki 



Applying this result to the expression of Eq. (44) we obtain 
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o 

where we have, in the last line, reduced the integration interval to half the 
original one. 

Finally, from the previous expressions for the axially symmetric case, the 
B(E2, liOi — > I/O"/) in the framework of the AMPGCM can be written as 



B{E2 , Uoi -> Ifa f ) = 



e 2 
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